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FOREWORD 


The present report contains the results of an investigation of ESSA 
VII Satellite radiation data for use in long-term earth energy experiments, 
performed under Langley Research Contract No. NAS1-11871 for the National 
Aeronautics and Space Administration. 

This research study, performed by Drexel University, is one part of 
a much larger effort by several institutions, including Colorado State 
University, Virginia Polytechnic Institute and State University, Link 
Temco Vought as well as cognizant personnel at NASA Langley Research Center. 
This team is studying satellite systems for performing long-term earth radi- 
ation balance measurements over geographical areas, hemispheres^ and the 
entire earth for periods of 10 to 30 years. A major portion of the total 
effort is responsive to the AA37E proposal and the proposed LZEEBE system 
(House and Sweet, 1973) . 

This investigation was conducted to provide necessary design data 
for use in developing the LZEEBE system. ESSA VII satellite employed 
plate and cone radiometers to measure earth albedo and emitted radiation. 
Each instrument had a black and white radiometer which discriTnl.nated the 
components of albedo and emitted radiation, and is identical in most 
respects to the proposed black and white balloons of the LZEEBE system. 

Earth measurements were made continuously from ESSA VII for ten months. 

The express purpose of this investigation is to obtain the ESSA VII 
raw data and process it to a point where it can be further analyzed for: 

• development of long-term earth energy experiments 

• document climate trends 

In addition to satisfying the immediate requirements of the contrac- 
tual statement of work, the research effort has included other areas of 
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FOREWORD (Continued) 


investigation. This report is one of three companion reports which together 
constitute the final report of phases X and II of subject contract. The two 
other reports are entitled; "Steady-State Solution to the Conduction Problem 
of a Spherical Balloon Radiometer," published as NASA CR.—132624, and "Our 
Contaminated Atmosphere - The Danger of Climate Change," published as NASA 
CR-132625. 

Gratitude is extended to several NASA/LaRC personnel for their en- 
couragement, interest, stimulating discussins and suggestions provided 
during the present investigation. Among these scientific personnel are 
included: Messrs. George Sweet (technical monitor), Charles Woemer 

and Jack Cooper, 
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INTRODUCTION 


Radiation balance measurements from satellites are worthwhile scien- 
tific experiments to perform. The radiation balance, defined, as the dif- 
ference between absorbed solar radiation and emitted thermal radiation at 
the outer boundary of the earth-atmosphere system, is the source of power 
that drives the atmospheric heat engine. Temporal and geographical dis- 
tributions of this net heating and cooling produce a broad spectrum 
of weather events. For example, long-term changes in the global radia- 
tion balance result in climatic variations. Heating in equatorial re- 
gions and cooling in polar regions cause thermal contrasts at the surface 
and in the atmosphere which produce seasonal changes. The radiation 
balance is one diabatic component of atmospheric motion which influences 
general circulation patterns over time scales of several weeks. Geogra- 
phical distributions of the radiation balance may affect the development 
of smaller scale features such as synoptic and mesoscale weather events. 
Observations of the global distribution of this differential heating and 
cooling at thr outer boundary of the atmosphere represent valuable scien- 
tific information. 

The net rate of radiant power gain for an area element of the earth 
on a specified time scale is determined by the relation: 

On ” H s - W r + V " (1 “ A > H s “ W 

where the parameters are shown schematically in Figure 1. All symbols 
in equation (1) and others to follow are defined in the List of Symbols 
in the front matter of this report. 

When integrated over the earth, describes the radiant energy ex- 
change between the planet earth and space. This quantity is determined 
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EARTH IRRADXANCE RADIANT EXITANCE 

INSOLATION REFLECTED EMITTED 



Figure 1. Radiant power balance on an earth surface element. 
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by the fluxes of incoming solar electromagnetic radiation, H , reflected 

solar radiation W , and emitted long-wave radiation, W , crossing a 

horizontal element of area just outside the earth's atmosphere (P = 0), 

2 2 

Representative values of these fluxes are*. 1353 w/m for H Q , 240w/m 

2 

for W , and about 400w/m for W at a location where <f> - 0. Other 
e r 

possible radiant power sources such as the fluxes of particle radiation 
and of electromagnetic radiation from the moon and all other stars, and 
the heat flux from the earth's interior are comparatively small. 

Satellite experiments are designed to measure the following; 

* The infrared radiant exitance, W , and 

C 

* The planetary albedo, A; where A = W r /H s , 

Satellites offer an ideal platform from which to perform radiation 
balance measurements owing to their locations outside the outer boundary 
of the atmosphere. In addition, polar orbiting satellites provide mea- 
surement coverage of the entire earth, but at nearly the same local time 
during the day or night. 

Notwithstanding the above attractive possibilities of performing 
radiation balance measurements from satellites, certain desirable appli- 
cations of the data are somewhat elusive. Problems such as the frequency 
and geographical distribution of sampling, spatial resolution of the mea- 
surements, and experimental accuracy pertain to both the design of experi- 
ments and the scientific applications of the data. As Moller and Raschke 
(1969) have pointed out, high accuracies of the components in equation (1) 

may be required. For example, if the global mean net radiation is 0.7 
2 

watt/m and an accuracy of 20 percent is desired, the global means of 

and (H - W ) are required to within an accuracy of 0.3 percent, 
s r 
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Obtaining these accuracies and properly sampling the earth both tem- 
porally and with sufficient spatial resolution are exceedingly diffi- 
cult, even with highly sophisticated instrumentation and a multiple 
satellite system. 

The foregoing discussions involving the applicability and associ- 
ated problems of observations from satellites has been presented so that 
the radiation balance measurements from ESSA weather satellites can be 
viewed in proper perspective. Instrumentation on these satellites pro- 
vide observations of W , W and an estimate of the difference (H - W ) . 

r e s r 

The radiometers have broad-band spectral and low-resolution spatial 
characteristics. The parametric measurements of radiant power density 
at satellite altitude involve an integrated average of irradiances from 
the earth-viewed area, as influenced by the optical characteristics of 
the radiometers and the rotation of the satellite. These measurements 
can be scaled to any arbitrary reference spheroid to account for varia- 
tions in satellite altitude. They are useful in the study of the global 
radiation balance and can probably be employed effectively in studies 
related to general circulation patterns. 

The major objective of the current research effort is to obtain 
ESSA VII raw data and process it to a point where it can be further 
analyzed for the development of long-term earth energy experiments and 
documentation of climate trends. Specific tasks in this regard are out- 
lined in the following statement-of-work; 

1. Obtain ESSA VII tapes from NESS/NOAA, 

2. Develop computer software to preprocess the data tapes. 

3. Determine inflight calibration constants, correlate them with 
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preflight calibrations and apply the results to the raw data, 

4. Perform color corrections to the white detectors, 

5. Correlate plate and cone measurements, 

6. Perform an error analysis. 

7. Develop a statistical technique for performing co-spectrum and 
cross-spectrum analyses, 

8. Insure that magnetic tapes can be processed with LRC computers. 

9. Evaluate the accuracy of inflight methods for calibrating 
radiometers and color constants, 

10, Provide guidelines for system design and flight operations of a 

proposed satellite system for measuring the long-term earth energy 
budget. 

It should be emphasised that the procedures discussed in the follow- 
ing sections of the report find direct application to the balloon radio- 
meters of the LZEEBE system. In fact, the geometrical viewing considera- 
tions of a spinning flat plate on ESSA VII are more complicated than the 
omni-directional spherical geometry of LZEEBE. 

A description of ESSA VII radiometers is presented first before 
discussing processing aspects of ESSA VII data. This is followed by a 
description of preflight instrument calibration and then a new method 
for inflight calibration of shortwave color constants, developed ex- 
pressly with the LZEEBE system in mind* using ESSA VII data. Comments 
concerning the correlation of cone and plate measurements follow the 
above sections. Preliminary conclusions concerning system design and 
flight operations of measurement systems conclude the main text of the 
report. 
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ESSA VII INSTRUMENTATION 


Historically, the measurement concept employed on LZEEBE and ESSA 
VII satellites is similar to earlier instrument packages f lown on Explorer 
VII satellite and TIROS III, IV and VII satellites (Suomi, 1958; Suomi, 
1961) . Two radiometers having different optical surface properties are 
used to discriminate the upwelling earth flux into its shortwave, W and 
longwave, components. Exposure of both radiometers to direct solar 
irradiance provides additional information concerning the magnitude of 
(H - W ) in equation (1) . Employing various combinations of these mea- 
surements, it is also possible to determine the earth’s albedo and the 
radiation balance parameter, providing assumptions concerning the angular 
variation of earth radiant exitance are made. 

Measurement theory for calculating the above parameters from tele- 
metry data has been discussed by several authors (Suomi, 1961; Weinstein 
and Suomi, 1961; Bignell, 1962; House, 1965; MacDonald, 1966; Suomi, et. 
al., 1967; Vonder Haar, 1968; House, 1968; House, 1970). These mathema- 
tical techniques apply to either mirror-backed hemispherical radiometers, 
flux plate radiometers (both mounted on a rotating satellite frame of 
reference) and isolated, omni-directional spheres of the LZEEBE system. 

The following discussion of this section includes the radiometer 
construction and viewing characteristics of ESSA VII instrumentation, 
and the general equation of absorbed irradiance on a unit detector area. 
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RADIOMETER CONSTRUCTION AND VIEWING CHARACTERISTICS 

A detailed description of the instrument package on ESSA weather 
satellites has been presented in detail by Nelson and Parent (1967) . 
Briefly, the experimental configuration involves a set of four radio- 
meters on each of two opposite sides of the satellite where each de- 
tector is connected electrically in series with its mate on the other 
side. The detector arrays are mounted on the circumference of the 
cartwheel satellite, perpendicular to the spin axis which is normal to 
the orbital plane (see Figure 2) . The temperatures of the detectors 
and housings supporting the detectors are the measured signals at the 
spacecraft. 

Each detector alternately views the earth and space during a rota- 
tion. Since the rotational period is relatively slow, about once every 
6.5 seconds, a single detector would have a modulating signal of about 
1°C to 3°C in phase with the rotational period. Since the integrating 
time of a temperature observation is relatively short compared to the 
rotational period, the thermal variation magnitude is significant in 
terms of earth flux calculations. However, the design of the experi- 
ment is such that the matched detector on the opposite side compensates 
for the modulating signal. As a result, the measured signal (temperature) 
is treated as being independent of the- spacecraft's rotation and is 
interpreted as an average temperature over the time span of the rota- 
tional period. 

Each of the two opposing arrays contain four flux plate detectors 
which are thermally isolated from the spacecraft as effectively aS possi- 
ble. Two of the four detectors in an array are exposed to a hemisphere 
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Figure 2. Conceptual drawing of ESSA VII satellite showing cart- 
wheel configuration, orbital plane and earth-viewed area. 


field-of-view (Figure 3a) whereas the other two detectors are located 
at the apex of highly reflecting cones (Figure 3b). For convenience, 
these radiometers are referred to as "plate 11 and "cone" radiometers, 
respectively. 

Each detector pair has different optical surface properties in order 
to discriminate the radiation balance components. One surface is painted 
black which has the property of being equally sensitive to both shortwave 
and longwave radiation. The other detector of the pair is an anodized 
aluminum or "white" surface which has the property of reflecting short- 
wave radiation and absorbing longwave radiation. Both of the cone and 
plate radiometers have black and white optical surfaces. 

ABSORBED IRRADIANCE ON A UNIT DETECTOR AREA 

During the course of a satellite orbit, between one and three sources 
of irradiance may be absorbed by the radiometers. If the satellite is iti the 
earth’s shadow (nighttime mode) only one source is present, W^. Just before 
ingression into, and egression from, the earth’s umbra (calibration mode) , two 
sources are present, W and H . During the remainder of the orbit (day- 

6 S 

time mode) all sources are present, H s and VT.. 

Since the instrumentation employed on ESSA VII comprise radiative 
equilibrium bolometers, all absorbed irradiances during the satellite or- 
bit are balanced by an equal radiant power loss per unit detector area 
The general mathematical expression for this balance is 

5 F W + a F H + a F W - N . ' (2) 

eeesssrrrd 

The three terms to the left of equation (2) comprise longwave emitted, 
direct solar and reflected shortwave irradiances, respectively. 


9 




Figure 3. Expanded view of ESSA VII plate radiometer components 
the radiometer with cone optics '(b). • Taken from Heist 
Parent (1967) . 


'a) and 
i and 
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For the plate radiometers on ESSA VII, the first term applies during 
nighttime mode, the first two terms apply during calibration mode near 
satellite ingression and egression, and all three terms apply during 
daytime mode observations. On the other hand, the optics of the cone 
radiometers shield the detectors from direct solar irradiance. Thus, 
only the first and third terms of equation (2) apply during the day- 
time mode. 

It should be clearly emphasized that some rather drastic assump- 
tions are being made when equation (2) is applied to real radiation ba- 
lance measurements. It is assumed that the optical properties of the 
detector over both the broad shortwave and longwave regions of the 
spectrum can be approximated by average absorptivities or and a^, 
respectively. For this to hold, the optical surfaces must have uniform 
spectral response over both broad-band regions since spectral irradiances 
in both regions vary as a function of wavelength. If this is not the 
case, some error will be introduced into the calculated radiant power 

densities W and W . 

e r 

A second major assumption is incorporated in the use of the shape 
factor F. The shape factor is a geometric term coupling all the angular 
absorbing properties of the detector configuration as well as assuming 
knowledge of any angular characteristics of the earth radiation field 
impinging on the detector. Certainly, the shape factor will vary both 
in time and geographical location, especially F^ which incorporates the 
anisotropic reflecting properties of the earth And atmosphere. The 
validity of both assumptions concerning constant absorptivities and 
shape factors should be investigated in depth in future research efforts 



Returning again to equation (2) , the radiant power loss by a unit 
detector area is expressed as 

A ■ A A 

— = oT + M(T)T + R(oT - aT ) 
e c 

T is the measured temperature of the radiometer and is the case tem- 
perature supporting the radiometer. These two temperatures are tele- 
metered to earth. The mass constant M(T) and the case radiation con- 
stant R are subjects of inflight and preflight calibration procedures, 
respectively, to be discussed in sections to follow. 

PROCESSING ESSA VII RAN DATA 

One of the major accomplishments during this investigation was the 
successful processing of 20 reels of ESSA VII raw data. It required 
about nine months of effort to overcome many unexpected pecularities in 
the data format. Certainly the documented format was far different than 
that which was really on the reels of tape, especially the positioning 
of end-of-file indicators. Many times these would occur in the middle of 
a satellite readout or there would be two or three end-of-file indicators 
at the end of each readout. One cannot appreciate the exasperation of 
processing data until he has attempted it himself. 

Part of the data processing effort included the development of tables 
which convert raw counts from the satellite to appropriate calibrated 
temperatures. In order to accomplish this task, hundreds of calibration 
points were plotted for each radiometer and case. Smooth curves were 
drawn through the average trend of these points. In some cases, there 
were obvious errors in the calibration data. All such points which were 



completely inconsistent with the other data were discarded. 

Details of the software used in processing the ESSA VII data are 
presented in Appendix A. The processed data are stored on 20 computer 
tapes in FORTRAN readable format which can be utilized on any computing 
system. Copies of these data are available from Drexel University upon 
request. 

During the data processing, each satellite readout is analyzed for 
bad data and consistency with orbit ephemeris. The status of each read- 
out was listed, independent of writing information on the new tapes, so 
that a formal documentation of data quality is available for future pro- 
cessing, For example, if a readout was found to be inconsistent with the 
ephemeris, the processing of that readout was terminated and an appropii- 
ate note of termination was printed on the data quality documentation 
listing. Then processing of the next readout was started. 

A tabulation of ESSA VII data available for future processing is pre- 
sented in Table 1, Useful data are available during the period 3 Sept. 
1968 - 22 June 1969. Evidently there was a problem with the tape recorder 
during the latter stages of data acquisition. There were several mis- 
matches of observations and ephemeris information, starting on tape #1018. 
Tape #1019 has only 4 days of data. There is an eleven day gap in data 
between the end of #1019 and the start of tape #1020 which contains the 
.last eight days of data. 

The detailed information presented in Table II tabulates specific 
readouts that were judged not suitable for future processing and also 
accounts for gaps in the data. Each tape covers about a period of one 
half month. There were several bad readouts on the first tape #1001 and 
on tapes #1015 - #1020, The record shows that data recovery diminished 
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during the last 2 1/2 months of operation. However, the amount of use- 
able data for future investigations is quite good when one considers all 
the number of useable readouts. During the entire period of data, there 
is 88% useable data. Of the remaining 12%, about 2% comprise inconsis- 
tent orbits and 10% of the readouts are missing. If one considers the 
time period between 3 September 1968 and 31 Hay 1969 , the useable data 
comprises 93% of the period. It must be remembered that a small por- 
tion of each readout is lost during each interrogation from the ground 
control. 



TABLE 1 


Available ESSA VII Satellite Data 


Tape ID 

Calendar Dates 


Julian Days 

n V.n r? 

Orbit Numbers 
Benin End 

1001 

3 Sept, 

- 17 Sept. 

'68 

18 

32 

231 

400 

1002 

17 Sept. 

- 30 Sept. 

'68 

32 

45 

400 

570 

1003 

30 Sept. 

- 15 Oct. 

•68 

45 

60 

570 

758 

1004 

15 Oct. 

- 31 Oct. 

*68 

60 

76 

750 

954 

1005 

31 Oct. 

- 15 Nov. 

'68 

76 

91 

954 

1146 

1006 

15 Nov. 

- 30 Nov. 

'68 

91 

106 

1146 

1334 

1007 

30 Nov. 

- 15 Dec. 

'68 

106 

121 

1334 

1521 

1008 

15 Dec. 

- 31 Dec. 

*68 

121 

137 

1521 

1722 

1009 

31 Dec. 

- 15 Jan. 

*69 

137 

152 

1722 

1910 

1010 

15 Jan. 

- 31 Jan. 

'69 

152 

168 

1910 

2110 

1011 

31 Jan. 

- 15 Feb. 

'69 

168 

183 

2110 

2298 

1012 

15 Feb. 

- 28 Feb. 

*69 

183 

196 

2298 

2458 

1013 

1 Mar. 

- 15 Mar. 

T 69 

197 

211 

2462 

2645 

1014 

15 Mar, 

- 31 Mar. 

'69 

211 

227 

2645 

2849 

1015 

31 Mar. 

- 14 Apr. 

’69 

227 

241 

2849 

3029 

1016 

15 Apr. 

- 30 Apr. 

'69 

242 

257 

3035 

3222 

1017 

30 Apr. 

- 15 May 

*69 

257 

272 

3222 

3409 

1018 

15 May 

- 31 May 

’69 

272 

288 

3409 

3611 

1019 

31 May 

- 4 June 

’69 

288 

292 

3611 

3655 

1020 

15 June 

- 22 June 

*69 

303 

310 

3797 

3888 
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.TABLE 2 


Inconsistent or Missing ESSA VII Satellite Data 


Tape ID 

Julian Day 

Orbits 

Status of 

Data 

Number of 




Inconsistent 

Missing 

Lost Orbits 

1001 

20 

254-258 

X 


2 


22 

270-274 

X 


4 


22 

278-282 

X 


4 


22 

282-287 

X 


5 


26 

328-332 

X 


4 


28 

353-357 

X 


5 


29 

366-370 

X 


4 


30 

370-374 

X 


4 


31 

384-388 

X 


1 


31 

388-391 


X 

3 

1002 

37 

462-466 

X 


4 


43 

537-541 


X 

4 

1003 

53 

668-674 


X 

6 


54 

684-689 


X 

5 

1004 

62 

779-783 


X 

4 


64 

804-813 


X 

9 


65 

817-821 


X 

4 


75 

946-950 


X 

4 

1005 

79 

988-991 


X 

3 


87 

1096-1102 


X 

6 

1006 

100 

1255-1259 


X 

4 


104 

1304-1309 


X 

5 

1007 

106 

1334-1338 


X 

4 


110 

1376-1380 





120 

1509-1513 




1008 

125 

1571-1576 

- 

X 

5 


126 

1580-1582 


X 

2 


. 136 

1705-1709 


X 

4 

1009 

143 

1798-1801 


X 

3 


145 

1826-1830 


X 

4 


16 



TABLE 2 - continued 


Tape ID 

Julian Day 

Orbits 

Status of Data 

Number of 




Inconsistent Missing 

Lost Orbits 

1010 

157 

1973-1977 

X 

4 


168 

2102-2105 

X 

3 

1011 

173 

2173-2175 

X 

2 


182 

2285-2289 

X 

4 

1012 

189 

2371-2374 

X 

3 


192 

2410-2415 

X 

5 

1013 

199 

2503-2507 

X 

4 

1014 

218 

2731-2737 

X 

6 


223 

2798-2801 

X 

3 

1015 

230 

2879-2884 

X 

5 


231 

2891-2897 

X 

6 


234 

2941-2947 

X 

6 


236 

2960-2962 

X 

2 


238 

2979-2991 

X 

12 


239 

3004-3006 

X 

2 


241 

3029-3035 

X 

6 

1016 

242 

3041-3046 

X 

5 


244 

3061-3063 

X 

2 


246 

3080-3085 

X 

5 


248 

3104-3110 

X 

6 


249 

3117-3122 

X 

5 


249 

3126-3129 

X 

3 


253^ 

3167-3172 

X 

5 


253 

3176-3179 

X 

3 


254 

3188-3197 

X 

9 

1017 

257 

3225-3228 

X 

3 


258 

3238-3242 

X 

4 


261 

3267-3272 

X 

5 


266 

3335-3338 

X 

3 


268 

3367-3373 

X 

6 

1018 

274 

3439-3443 

X 

4 


279 

3489-3492 

X 

3 
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TABLE 2 - continued 


Tape ID 

Julian Day 

Orbits 

Status of 

Data 

Number of 




Inconsistent 

Missing 

Lost Orbits 


279 

3498-3504 


X 

6 


280 

3504-3505 

X 


1 


283 

3542-3548 


X 

6 


285 

3573-3576 


X 

3 


287 

3598-3601 

X 


3 


288 

3605-3611 

X 


6 

1019 

292 

3655-3797 


X 

142 

1020 

304 

3806-3810 

X 


4 


307 

3849-3851 

X 


2 


309 

3880-3882 


X 

2 


Summary of Data Recovery 

First Orbit - #231 

Last Orbit - #3388 

Total Possible Number - 3657 Orbits 

61 Inconsistent Orbits 

376 Missing Orbits 
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PREFLIGHT INSTRUMENT CALIBRATION 

Several constants must be known in equations (2) and (3) before 
calculations of earth radiation parameters can be accomplished. These 
include the average broad-band absorptivities in equation (2), and long- 
wave emissivity, mass and case radiation constants in equation (3) . All 
of these terms in the equations should be determined prior to launch and 
then checked thru inflight calibration procedures. However, in the case 
of the ESSA VII radiometers, the only preflight laboratory calibration 
which took place was enough data to determine the case radiation constants 
for the cone and plate radiometers. In addition to this information, data 
concerning the longwave optical properties of the black and white radio- 
meters can be used to establish a and e. 

Li 

Results of the investigation concerning the absorbing and emitting 
properties of the ESSA radiometers as well as the computation of the case 
radiation constants are presented in this section of the report. 

LONGWAVE COLOR CONSTANTS 

As a first approximation, the longwave absorptivities and emissivities 
are computed from theoretical models using available spectral absorptance 
measurements of the sensors' surfaces. It should be emphasized that these 
measurements may be somewhat different in absolute magnitude than the ab- 
sorptance properties of the flight models. However, It is reasonable to 
assume that the gross variation of the absorptances with wavelength is 
similar; it is this variation that is being evaluated here. 

Equations for computing absorptivities and emissivities are as 
follows: 
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and 


- _ 


50u 

j a(X) N(X) dX 


a = 


50u 

l 

Ijj 


N(X) dX 


L 


50u 

e (X) B(T, X) dX 


ly 

50y 


L 


B(T, X) dX 


( 4 ) 


(5) 


where from Kirchoff's Law, e(X) a a(X) , It is implicitly assumed here in 
the calculations that the observed spectral absorptances of each surface 
are valid for all temperatures of the detectors. This, In fact, may not 
be a valid assumption and should be investigated further. 

When computing absorptances using equation (4) , it is necessary to 
specify the variation of the earth flux with wavelength. A simple model 
is assumed which states that the variation In N occurs only in the atmos- 
pheric window region, between 8 and 13 microns. All other wavelengths 
emit at a constant temperature of 216. 6°K which is the tropopause tempera- 
ture of the Standard Atmosphere. This model grossly simulates the spectral 
emitting properties of the earth-atmosphere system. However, the results 
here agree with the correlation of nighttime plate observations presented 

later. In mathematical terms, this model states: 

- j B (216.6°K,X) for X < 8p and X > 13u 
|B(T,X) for 8p - X ^ 13n 


20 


Calculations of the absorptivities and emissivities were performed 
for a variety of temperatures and simulated earth fluxes using equations 
(4) thru (6), and using the spectral absorptance data in Figure 4„ The 
results of these calculations are presented in Figures 5 and 6 for the 
black and white radiometers, respectively. These results are plotted in 
terms of equivalent black body temperature which would correspond to the 
detector temperature in the case of emissivity, and the equivalent black 
body temperature for the longwave radiant exitance for the absorptivities. 
It is evident from these results that broad-band optical properties 
may vary significantly with temperature, even for the black paint. The 
uncertainties in this conclusion are the unknown variation of the spec- 
tral absorptances with detector temperature and the true character of the 
spectral absorptance curve beyond a wavelength of 25 microns. One cer- 
tainty here is that proper calibration of the spectral properties of opti- 
cal surfaces are important to the success of any future instrumentation 
such as the LZEEBE satellite system . 

£n order to process ESSA VII radiation data, some selection of long- 
wave optical properties must be made from the information contained in 
Figures 5 and 6. In this regard, it is reasonable to select color con- 
stants for the nighttime and daytime mode of calculations. These selected 
values are presented in Table 3. It will be shown later that the selection 
of longwave color constant is consistent with nighttime observations (see 
section concerning correlation of plate and cone measurements) . 

POWER EXCHANGE BETWEEN RADIOMETER AND CASE 

Preflight laboratory calibration data were available which provide 
information for determinin g case radiation constants for all radiometers 
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M^UTZ 9300-S FLAT BLK PAINT (Schoffer, 1962] 


{-—Extrapolation to 50/i 



ANODIZED ALUMINUM 
SAMPLE #*7 (Sparkman, 1964) 


-H K- AVERAGING INTERVAL \p. 


WAVELENGTH - MICRONS 


Figure 4. Spectral absorptance properties of the black artd white ESSA 
radiometers. 





MAUTZ 9300 -S FLAT BLACK PAINT 


ABSORPTIVITY 


. EM1SS1V1TY 


EQUIVALENT BLACK BODY TEMPERATURE, °K 


Figure 5 . Variation in absorptivity and eraissivity with temperature 
black' radiometers . 







EQUIVALENT BLACK BODY TEMPERATE 


Figure 6. Variation in absorptivity and emissivity with temperature 
white radiometers . 





TABLE 3 


Longwave Color Constants for ESSA VII Radiometers 


Radiometer 

ae 

Hade 


n 

Nighttime e^ 

Daytime 


White Plate 

.80 

.74 

.76 

0.97 

Black Plate 

.94 

.92 

.94 

0.98 

White Cone 

.80 

. .75 

.75 

1.00 

Black Cone 

.94 

.93 

.93 

1.00 




on ESSA VII. The calibration procedure consists of placing the flight 
version of the sensors in an enclosed vacuum chamber whose walls are main- 
tained at a constant temperature. For all practical purposes, the chamber's 
surface is assumed to emit as a blackbody since its surface is painted 
black, and its area is large compared to the area of the sensor's surface. 
The chamber is maintained at a cold temperature (around -65°C) and its 
heat capacity is quite large compared to that of the sensor. Thus, the 
chamber maintains a known radiation field within the confines of the en- 
closure . 

The sensor's case is heated electrically, and both the sensor and 
mount temperatures are monitored. These temperature data and the chamber 
temperature provide the available information for calibrating the sensor 
system. 

The case radiation constant can be computed from an expression de- 
veloped from equations (2) and (3) . Since the radiant field within the 

enclosure is assumed to be black body * radiation, the irradiance on the 

4 

radiometers is oT where T is the temperature of the source (i.e., the 
s s 

chamber) . In radiative equilibrium, power lost by the radiometer back to 
the chamber is W^. Equating these two components gives the expression 


p 4 4 4 4. 

— oT = oT + R(oT - a T ) 
e s c 

n 


Solving for the case radiation constant yields 


n 



< 7 > 


CS) 
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Incorporating longwave color constants listed in Table 3, calculations 
of case constants were performed using equation (8) . The results of 
these calculations are presented in Table 4. 


Case Radiation 

Radiometer 
White Plate 
Black Plate 
White Gone 
Black Cone 


TABLE 4 


Constants for ESSA VII Radiometers 


Case Constant R (no units) 
0.075 
0.054 
0.36 
0.33 
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INFLIGHT CALIBRATION PROCEDURES 


The previous section of the report discussed the determination of 
the longwave absorptivities and emissivities , and the case radiation 
constants for ESSA VII radiometers. Two parameters remain to be deter- 
mined - the mass constants and the shortwave color constants. These 
parameters are best determined in the space laboratory where there is a 
hard vacuum and the best possible source of shortwave irradianee, the sun. 
The observational data themselves are used to determine the magnitudes of 
these parameters, utilizing the sun as a calibrating source of shortwave 
radiation. 

The accurate determination of the shortwave color constants and the 
accurate monitoring of any changes in their values are of utmost impor- 
tance to the confidence in ESSA VXI data as well as the observational 
integrity of the LZEEBE system . This statement cannot be emphasized too 
many times. Therefore, every effort should be exerted to perfect the 
inflight calibration procedures. 

Data utilized for these purposes are acquired during thermal tran- 
sitions at ingression cool-down and egression warm-up of the radiometers 
exposed to direct solar radiation. At these particular locations of the 
satellite orbit, reflected earth radiation can be neglected. Thus, 

the only sources of irradiation are H and W . If accurate calculations 
J s e 

of W can be determined at times just before ingression and/or just after 
egression, it is possible to determine the radiometers’ response to short- 
wave irradiance by subtracting the W magnitude from the observations. 
Thus, the radiometers can be calibrated, using the sun as the source, if 
the longwave earth term can be determined. 
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In order to perform the above shortwave calibration procedure, it 
is essential that the mass constants of the radiometers must be accurately 
known so that magnitudes of can be calculated during thermal transi- 
tions. The goals then of inflight calibration procedures are to accom- 
plish three things! 

f accurate measurement of during thermal transitions. 

* frequent calibration of shortwave color constants. 

* periodic calibration of mass constants. 

A considerable amount of effort was directed to this end, using ESSA 
VII data. Needless to say, several rather sophisticated approaches were 
tried but did not return good results. In the end, a simple approach 
turned out to be the best procedure. This approach is presented below. 

CALIBRATION OF MASS CONSTANTS 

Consider first the radiative balance equation just before ingression 
or after egression, i.e., calibration mode, 


A # A 4 v 

— - F W + — =■ F H = cT + M(T)T + R(oT - aT j 
e e e e s s c 

When the spacecraft passes into the earth's umbra, H g is aero, and the 

equation takes on the form typical of nighttime calculations. 

— - F W = OT^ + M(T)T + R(crT^ - aT 4 ) 
e e e c 

If the terms in equation (10) are arranged in the following manner where 

4 4 , 4 4 \ 

oTj = aT + R(oT - aT ) 

1 c 

, 5 

aT? = M(T) (-T) + — - W 

i e ee 

This equation has the form of a straight line, y - + Cq. 


(9) 


( 10 ) 


( 11 ) 
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In this regard, the slope of the straight line is the mass constant 
and the intercept Cq is the longwave earth term. This suggests that 

4 

thermal transition data should be plotted on linear graph paper as oT^ 
vs . (-T) . By analyzing the characteristics of this plotted function it 
is possible to establish the magnitude of the mass constant for the ra~ 
diometer. Furthermore, changes in the slope of the function would sug- 
gest a variable mass constant. 

Unfortunately life is not this easy. Implicit in the above discus- 

a 

sion is the assumption that the longwave term — did not change 
during the thermal cool-down. This would be true if the satellite were 
stationary for about 6 minutes during coolr-down. In real operation, low 
altitude spacecraft move about 1300 km ground track during this transi- 
tion. As a general rule the earth term will not he constant during the ac- 
quisition of cool-down data. As a result, this will lead to distortions 
of the straight line equation discussed in equation (11) . If a straight 
line were fit to the observed data through regression techniques, the 
functional distortions would manifest themselves in terms of slope changes. 

In order to obtain a better understanding of the problem Luther 
(private communication) provided ingression cool-down data for a simu- 
lated LZEEBE balloon. The data are idealized in that the balloon has a 
uniform temperature at any time during thermal transition; however, his 
data are si mil ar to measurements received from ESSA VII satellite. Three 

cases were analyzed with different longwave gradients during cool-down. 

* 

In the control case, W was constant, i.e., W =0. In the other two 

' e e 

^ 2 
cases, = + 0,2 (watts/m )/sec. 
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The results of the analysis of simulated balloon data are presented 

in Figure 7 where the normalized radiant power is plotted as a function of 

(-T). It Is evident from this plot that the data where the longwave earth 

term was constant appears as a straight line. When there is a significant 

gradient in the earth term the function is concave downward for W g <0, and 

concave upward for W >0. In other words, longwave gradients during cool— 

e 

down distorted the linearity of the function. In addition, the average 
slopes of the functions increased for W <0, and decreased for W >0. 
Analysis of simulated data such as this was a tremendous help in under- 
standing observations from ESSA VII satellite. 

Based upon these results, it was decided to filter the cool-down data 
from ESSA VII and select only those cool-down functions which exhibit li- 
nearity. In fact, such a linearity criterion could he incorporated in the 
electronics logic aboard the spacecraft in order to select useable cali- 
bration data. These data selections were further verified by examining 
the cone temperature data during the same temperature transitions. 

A similar procedure as outlined above can be followed in analyzing 
warm-up data during satellite egression. For these data, the variable 
term affecting functional linearity is still associated with longwave 
earth radiation. There is also a large constant associated with direct 
solar irradiation (see equation 9) , which results in a large value for the 
intercept value Cq in a curve-fit of the data. The slope of the line in 
a linear curve- fit of the egression is still related to the mass constant, 
A selection of data was taken from filtered ingression and egression 
cases throughout the 10 month period of ESSA VII measurements. 


31 



Each of these cases was fit with a straight line, and then the data 
were plotted on one graph in a normalized fashion as shown in Figure 7. 
Ingression and egression data were treated separately. 

A detailed study of the data showed that a smooth curve drawn 
through the mass of points still exhibited a change in slope during 

both the cool-down and warm-up transitions. This change in slope is 

* 

attributed to a variable mass constant. A cross plot of T vs, T was 
generated from the data in order to relate the variable mass constant 
in terms of radiometer temperature. The results of all this effort are 
presented in Figure 8 where the mass constants are plotted as functions 
of plate temperatures. As can be seen from the graphs, the mass con- 
stants are no longer constant, but linear variables increasing as a func- 
tion of radiometer temperature. A more appropriate name for H(T) might 
be the "mass term." 

The uncertainties of each of the points are indicated in Figure 8. 
The vertical uncertainty in mass term magnitude is based on the ability 
to determine the slopes of the ingression - egression data discussed in 
the previous paragraph. The horizontal uncertainty in temperature is 
based on the scatter of points in the cross plot of T vs. T. 

It should be mentioned here in a historical sense that there were 
several fruitless attempts to determine any significant variations in the 
mass constant as a function of temperature. The key to the solution of 
the problem was data filtering, i.e. , selecting only those cases that 
exhibit linearity and thus some assurance that the longwave earth radia- 
tion term was constant during thermal transitions. 
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© INGRESS !0N 
B EGRESSION 



CALIBRATION OF SHORTWAVE COLOR CONSTANTS 


The problem of calibrating shortwave color constants rests primarily 

a 

e 

on the ability to calculate accurate longwave earth terms — - — F^ W 
during the thermal transition. Since accurate values of the mass term 
were established in Figure 8, it follows that accurate values of the 
shortwave color constants can be calculated. 

During the course of this portion of the investigation, numerous, 
sophisticated approaches were attempted with little results. After a 
careful examination of the failures of other approaches, the simplest of 
procedures seemed to work best. 

The concept of the procedure is presented in Figure 9. The solid 
points are the radiant power densities associated with an ingression 
cool-down of the black plate on ESSA VII. The assumed time of ingres- 
sion is indicated on the graph. There is a two minute time interval 
following this time which is a period of solar contamination, i.e., the 
sun's disc is setting beneath the viewed horizon of the earth at satellite 
altitude. This time period is associated with the penumbra of the earth’s 
shadow. Before this period daytime mode calculations are appropriate (see 
equation 9) - after this period, nighttime mode calculations are used (see 
equation 10) . The middle time of this 2 minute period is defined as the 
time of calibration (see t at the bottom of Figure 9). 

The last four daytime mode calculations (first four solid points on 
graph) are used to extrapolate the trend of observations forward in time 
over a one minute period to the time of calibration. In a similar manner, 
the first four computed longwave earth terms (see first four open points 
at bottom of graph) are used to extrapolate the trend line of the longwave 
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Figure 9 , Concept of shortwave color constant calibration procedure* 
(see text for description) 
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term backwards in time to the time of calibration. The computed magni- 
tudes of these two trend line extrapolations at the precise time of cali- 
bration are subtracted from each other. The resulting magnitude of radi- 
ant power density constitutes the response of the black plate solely to 
direct solar irradiation. In the example shown in Figure 9, the follow- 
ing magnitudes arise from this procedure: 

2 

F W = 402.5 watts/m 
e e 

2 

F W = 59.3 watts /m 

e e 

2 

F H = 343.2 watts/m 
s s 

The solar term is indicated on the graph in Figure 9 at the time of cali- 
bration. 

The above procedure works even better for egression warm-up data . 
since the radiometer approaches an equilibrium temperature more rapidly 
than in the case of cool-down. 

Determinations of the solar term from consecutive cool-down, warm-up 

thermal transitions indicates differences of magnitudes of less than one 
2 

watt/m on the average. Good agreement such as this leads one to have 
confidence in the shortwave calibration procedure and in the determina- 
tion of mass constants. 

ANALYSIS OF ERRORS 

The question always arises as to the accuracy of the inflight calibra- 
tion procedures. This pertains to both the absolute accuracy of the cali- 


a a 

— F H 4- — 
s. s s e, 
d d 


e a 

( — ) ~ 
£ d E n 


Thus, the solar terra becomes 
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brations as well as determining the change in the solar terms of each 
radiometer. It should be pointed out that the absolute temperature ac- 
curacies of the observations are not considered in the evaluation of 
the calibration procedure since they are assumed to be correct. 

Without question, the accuracy of the inflight calibration proce- 
dure depends on the accurate determination of the mass term in Figure 8. 
Within the uncertainties of the points plotted on the graph, one can de- 
termine a radically different change in mass term with temperature sim- 
ply by changing the slope of the line, and still remain within the un- 
certainties indicated. However, such a change will produce inconsis- 
tencies in the solar terms between ingression and egression thermal 
transitions. Thus, it is possible to set the slope of the mass term 
change with temperature with some certainty, based solely on inflight 
data . 

Another aspect of the evaluation of the curves in Figure 8 which 
pertain to calibration accuracies is the translation of the curve up or 
down, a bias change, without affecting the slope of the curve,.. This 
adjustment truly affects the absolute accuracy of the solar term ob- 
servations within the uncertainty of the points on the graph. This 
uncertainty amounts to + 2 mass term units for the black plate and 4- 3 
units for the white plate. This uncertainty results in about a + 1 watt 
uncertainty in the solar terms of each of the plate radiometers. The 
only way to reduce this error is to consider hundreds of cases in order 
to drive down the random uncertainty of the mass term determinations . 
Such an approach is certainly feasible. 
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When one considers other errors related to the accurate determina- 
tion of shortwave and longwave earth fluxes such as the spatial sampling 
problem, the evaluation of shape factors for the anisotropic reflection 
problem, and the assumption of constant broad-band absorptivities , one 
quickly concludes that these problems are larger in magnitude than the 
errors associated with inflight calibration of shortwave color constants. 
It is estimated that these other problems are in the one to two percent 

error range, whereas, the inflight calibration uncertainty on the order 
2 

of one watt/m , or a few tenths of a percent error. 

It is concluded from the discussion in this section of the report 
and from the investigation of the inflight calibration problem that 
valid shortwave calibrations can be obtained solely from inflight data 
with a sampling frequency every 30 seconds . 
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CORRELATION 0? PLATE AND CONE ItEASUREVENT R 


Correlations of plate and corte measurements during nighttime mode 
operation give considerable insight as to the relative longwave optical 
characteristics of the radiometers. During the nighttime mode, the pair 
of plate radiometers, and the pair of cone radiometers view identical 
longwave irradiance fields, By imposing this boundary condition, it is 
possible to evaluate the relative a^/e ratios of the white and black opti- 
cal surfaces. 

In a similar manner, correlation of plate and cone measurements give 
information concerning the lelative geometrical viewing properties of the 
two radiometer types. It must be remembered in the latter case, however, 
that the cones absorb longwave irradiance from a restricted portion of 
the earth viewed area. Thus, the boundary condition that the cones and 
plates view the same irradiance field does not hold exactly, but is a rea- 
sonable approximation for the mean of many observations. 

Correlation of plate and cone measurements from ESSA VII were per- 
formed during this investigation for many observations throughout the ten 
month period or useful satellite data. Typical correlations for a selec- 
tion of these data are shown in Figures 10 and 11. 

The correlation of the white plate radiometer with the black plate 

radiometer is shown in Figure 10 where the white power density is plotted 

on the ordinate vs. black power density on the abscissa. By plotting the 

measurements in this manner, the slope of the linear regression line through 

the points is equal to the ratio of ( a e J / { ~*e } . The slope 

\ e / white / \ e / black 

of the line for the data in Figure 10 indicates a magnitude of (1.06). 
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Figure 10. Correlation of white plate vs. black plate radiometers 
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This result suggests that the (a e /e) ratio of the white plate is roughly 
6 percent larger than a similar ratio for the black plate. Referring to 
the values of longwave optical properties presented in Table 3, the com- 
puted ratio for white vs. black plate, using nighttime values, is 1,058 in 
magnitude, in good agreement with the correlation of nighttime observations 
presented in Figure 10. Thus, one has confidence in the longwave optical 
properties developed in the section on preflight instrument calibration, 
even though a crude model was employed for the upwelling spectral irra- 
diance from the earth. 

It is noted further from the correlation in Figure 10 that the stan- 

2 

dard error of the estimate is 0.4 watts/m , that is, about 68% of the points 

2 

fall within +0.4 watts /m of the line. In order to interpret this re- 
sult in terms of W earth flux values, a multiplier of approximately 4 is 
e 

2 

appropriate, giving a + 1.6 watts/m standard error. This magnitude of un- 
certainty is rather large, about 0.7% of the average W magnitude. In fact, 
this uncertainty is larger in magnitude than the errors in the solar terms 
of each radiometer discussed in the previous section. The reason for this 
uncertainty is the point to emphasize here. It is due to the fact that 
the spectral absorbing properties of the optical surfaces are not uniform 
across the longwave, broad-band spectrum . This point is easily understood 
by considering the observed spectral absorptance properties of the anodized 
aluminum "white" sample shown in Figure 4. 

As mentioned earlier in this section, the correlation of plate and cone 
measurements gives information concerning the relative geometrical viewing 
properties of the two radiometer types. In this regard, it is assumed that 
the optical properties of the black paint or anodized aluminum surfaces are 
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the same - certainly a reasonable assumption. An example of this correla- 

tion is presented in figure 11 where the black plate observations are 

plotted as a function of concurrent black cone radiant power densities, 

A linear regression line was computed for these data, in a manner similar 

to the results in figure 10. The slope of this line has a different in- 

terpretation than in the plate correlation. The slope is related to the 

ratio of shape factors for the plate and the cone radiometers, i.e., 

slope = (F ) . /(F ) . This ratio is 0.72 for the data shorn in 

£ piflts £ cone 

Figure 11. 

The results in Figure 11 have further engineering design implications 

concerning the LZEEBE system. The standard error of the estimate for the 

2 

data shown in Figure 11 is + 1,1 watts/m , When this uncertainty is in- 

2 

terpreted in terras of W , a variation of about + A. 4 watts/m is obtained. 

Similar uncertainties were obtained in a prior study by House (1968) when 

he correlated the change in plate observations vs. the corresponding change 

in cone observations, utilizing ESSA III data. This magnitude of uncertainty 

is greater than 1% of the W irradiance field. Therefore, it is concluded 

6 

that auxiliary, longwave gradient detectors would not have significant value 
in defining the longwave earth term during inflight calibration of shortwave 
color constants . In fact, the error in the plate - cone correlation is 
larger than the uncertainty in the inflight calibration procedures developed 
in the last section. 

STUDY CONCLUSIONS RELATED TO LZEEBE SYSTEM CONSIDERATIONS 

An attempt should be made here at the end of this report to summarize 
the conclusions of the investigation and how these conclusions relate to 
the engineering design and flight operation of the proposed LZEEBE system. 
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vs.* black cone radiometers. 
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It may be recalled that the major objective of this effort was to 
obtain ESSA VII raw data and process it to a point where it can be fur- 
ther analyzed for the development of long-term earth energy experiments 
and documentation of climate trends. It is felt that this objeccive was 
accomplished during the reporting period. Data recovery during the 10 
month observation period of ESSA VII was about 90% efficient. Certainly 
there is an abundance of useable observations for future investigations. 
Data are stored on 20 reels of computer tape in a FORTRAN readable format. 

A considerable amount of effort was spent in developing successful 
preflight and inflight calibration procedures for determining critical 
constants required for the future calculation of radiation balance para- 
meters. Even though these procedures pertain to ESSA VII radiometers 
they find direct application to the balloon radiometers of the LZEEBE 
system. In some ways the ESSA VII data reduction problem is more compli- 
cated than that for the LZEEBE system since a spinning flat plate has a 
more complicated geometry than an omni-directional sphere. 

Specific conclusions of this study and their relationship to engi- 
neering design and flight operations of the LZEEBE system are outlined 
below: 

ENGINEERING DESIGN CONSIDERATIONS - LZEEBE SYSTEM 
Optical Constraints 

detectors should have uniform spectral response to both the 
shortwave and longwave broad-band spectral regions. 

knowledge of the spectral absorbing properties of optical 
surfaces must be known between about 0.2 microns to about 
50 microns or longer wavelengths, if possible, 

the variation in spectral absorbing properties must be known 
for the expected ranges of radiometer temperatures and for 
the above spectral, region. 
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white optical surfaces like anodized aluminum and teflon 

• backed with aluminum have inferior longwave optical pro- 
perties compared to paint binders. 

white paints with binders similar to the black paints may 

• be preferred to the above materials provided the shortwave 
optical degradation Is not too severe. 

Geometrical Constraints 


knowledge of tie angular absorbing properties of the opti- 
, cal surfaces with wavelength is desirable to develop accu- 
rate geometric shape factors. 

optical surfaces should be exposed to direct solar irradiance 
, uniformly, if possible, so that any changes in optical proper- 
ties will be the same over the balloon’s surface area. 


Other Considerations 

crude longwave gradient detectors will not provide the needed 

• accuracies of earth radiation gradients during thermal transi- 
tions - however, these detectors may be useful in determining 
rough satellite orientation to solar radiation. 

all physical and optical properties of the balloon radiometers 

• should be known prior to launch - this enhances the value of 
inflight calibration results. 


FLIGHI OPERATIONS - LZEEBE SYSTEM 


Type of Observation or Calibration 

Data Sampling 
Time Interval 

Time Frequency 
of Performance 

W e and W r (daytime mode) 

30 to 120 sec 

30 to 120 sec 

W (nighttime mode) 

30 to 120 sec 

30 to 120 sec 

W (calibration mode) 
e 

30 sec 

30 sec 

Mass term calibration M(T) 

30 sec 

all orbits - analyze 
filtered orbits 

ct 

Calibration of solar term , s ^ . 

ingression & egreoaion r s s J 

30 sec 

all orbits - analyze 
filtered orbits 


Albedo calculation A 


Calibration F H 
s s 


/s \ /« \ 

tM - r* 1 

\a / V a /, 
L\e/a \e^b 


system independent 
(use W r averages) 


30 sec 


over time period _ 
of data analysis 

all orbits - cross check 
with solar term calibra- 
tions 


46 



REFERENCES 


Bignell, K, J. } 1962: Heat balance measurements from an earth satellite - 

An analysis of some possibilities. Roy. Meteor. Soc . , 231-244. 

House, F. B., 1965: The Radiation Balance of the Earth from a Satellite. 

Ph.D. Thesis, Department of Meteorology, The University of Wisconsin. 

, 1968: The reduction, analysis and interpretation of radiant 

flux measurernen' s from ESSA weather satellites. Final Report under 
ESSA/NESC Contract No. E-108-67(N) , GCA-TR-68-20-G , GCA Corporation, 
Bedford, Massachusetts. 

, 1970: The reduction, analysis and interpretation of radiation 

balance measurements from ESSA Weather Satellites. Final Report under 
ESSA/NESC Contract No. E-119-68(N), GCA-TR-69-17-G, GCA Corp., Bedford, 
Massachusetts . 

, G. Sweet, et.al., 1973: Long-term Zonal Earth Energy Budget 

Experiment (LZEEBE) . A proposal to AAFE, April. 

Luther, M. , 1974: Private Communication 

MacDonald, T. , 1966: Unpublished seminar presented in November at the 

Department of Meteorology, The University of Wisconsin, Madison, Wise. 

Moller, F. and E. Raschke, 1969: Problems of meteorological observations 

from satellites. Space Science Reviews , Vol, 9, No, 1, 90-148. 

Nelson, D. F. and Parent, R r , 1967: The prototype flat-plate radiometers 

for the ESSA III satellite. Annual Report on Grant No. WBG-27, Depart- 
ment of Meteorology, The University of Wisconsin, 

Suomi, V. E., 1958: The Radiation Balance of the Earth from a Satellite, 

Annals of the IGY , VI, 331-340. 

, 1961: The Thermal Radiation Balance Experiment on Board 

Explorer VII Satellite, NASA Tech. Note D, 608. 

, Hanson, K.K. , and Vonder Haar, T. H. , 1967: The theoretical 

basis for low-resolution radiometer measurements for a satellite. Annual 
Report on Grant No, WBG-27, Department of Meteorology, The University of 
Wisconsin. 

Vonder Haar, T. H, , 1968: Variations of the Earth's Radiation Budget. Ph.D. 

Thesis, Department of Meteorology, The University of Wisconsin. 

Weinstein, M. and Suomi, V., 1961: Analysis of satellite infrared radiation 

measurements on a synoptic scale. Monthly Weather Review , 89, 419-428. 


47 



APPENDIX A 


COMPUTER SOFTWARE FOR PRELIMINARY PROCESSING 
OF ESSA VII SATELLITE DATA 
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DESCRIPTION OF THE MAIN PROGRAM 

A two block approach is used in the processing of all data. The main 
program calls the input procedure and it is from there that all of the 
data is read in. Data is read from the CDC tapes one block (25 lines) 
at a time and assembled into the two arrays ORBIT and SENSOR. These two 
arrays are then passed back into the main program. There are always two 
blocks of data in storage at any one time. Differences in angles and 
count values are computed in order to perform the necessary counts to 
temperature conversion, and as the last lines of each block are processed 
some overlap into the next block is necessary. Upon completion of all 
work on a block, that block (consisting of the finished ORBIT and SENSOR 
arrays) is written onto tape and the next block of data is read in. This 
approach reduces the total amount of time needed since computations on 
each block are being performed at the same time that the previous block 
is being written. Parts A and B of the main program segment are there- 
fore identical and control is passed to either part depending upon which 
particular block is being processed. 

The first loop in PARTA (down to label PHILA) determines the day-night 
position of the satellite and adds 1024 counts onto the black plate day- 
time count values. Bad rows are skipped. The satellite position logic 
is as follows. The principle testing point was the zenith angle (column 5 
of ORBIT array). Experience showed that the day-night transition occurred 
within a window of zenith angle values ranging roughly from 120 to 128 de- 
grees, with the angle increasing in going from day to night and decreasing 
from night to day. If the zenith angle was less than 115 degrees, daytime 
was assumed, and 1024 counts were added onto the black plate values. All 
other lines with zenith angles greater than 115 degrees were subject to 


49 


further testing. Ensuing tests required 1) the sense of direction of 
the zenith angle (either increasing or decreasing) and 2) the gradient 
of the black plate count values. Actual processing indicated the need 
for one additional test because in some readouts bad data occurred in the 
transition zone and the change from day to night (or vice versa) was not 
detected. A third variable was computed, the absolute difference between 
the black and white plate count values. These two values agree quite well 
during the nighttime observations and they differ greatly in daylight. If 
this difference was less than 100 counts then night was assumed. Finally, 
if white plate counts ran greater than 350, then daytime was assumed. On 
testing any one good line of data, if the succeeding line contained a 
black plate count value of 9999 then the line being tested was replaced 
with 9999 since the black plate gradient and other variables could not be 
computed. This results in the loss of some good data but the generally 
good condition of the data makes this loss negligible. The program segment 
from the label PHILA to the label IOWA performs all of the above operations 
using the last lines of one block and the first lines of the succeeding 
block. The segment from the label IOWA to the label ELM works n the first 
column of SENSOR array which contains the reference resistor count as well 
as the mount count values. Again bad data lines are skipped and since the 
count values of these sensors do not fall below 400 any values below 400 
have 1024 added to them. The segment from label PBY to label ILL converts 
count values to temperatures for the black and white plates and cones. All 
count values are adjusted using the reference resistor count (all counts 
are multiplied by .994 for ESSA 7). Since the temperature tables had to be 
split into two parts of 1000 and 500 lines, count values had to be tested 
and converted to temperatures using the appropriate table. In order to 
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safeguard against a stoppage in program execution due to negative count 
values (count values are a parameter of the temperature table arrays) any 
count values exceeding the acceptable limits were automatically set to a 
value of 500. This would effectively place a temperature of 0.00 at that 
location . An inspection of the temperature tables shows that no tempera- 
ture of exactly 0.00 can exist in the data. Therefore, in performing any 
additional processing of this data a test for zero temperatures on all 
seven sensors should be made. Appearance of these zero values will indi- 
cate original data which were out of the temperature tables range. The 
program segment from label TEX to label NH. uses the above procedures to 
convert mount count values to temperatures. The remaining lines from 
label MM to label PARTB write the finished data onto tape, write end of 
file on tape increments the record counter, and reads in a new block of 
data via an input statement. PARTB comprises the rest of the program and 
is identical to PARTA. 

INPUT PROCEDURE DESCRIPTION 

The input procedure is used to read in new data from the input tapes, 
screen the data, and set up the ORBIT and SENSOR arrays which are passed 
back into the main program. The twenty eight work documentation headings 
as well as most of the. end of files are also written onto the output tape 
from this procedure. When a readout is finished and an end of file written, 
the record counter in the main program is set to zero and the input proce- 
dure is called again. There are two read statements in this procedure 
labelled MINN and ONT respectively. MINN is executed when the record 
counter is zero and it is used to read in the initial input files consis- 
ting of the twelve word label as well as the twenty eight word documenta- 
tion headers . Upon writing the twenty eight word heading control Is 
passed back to the main program, the record counter is incremented, and the 
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input procedure is called again. Now data will be read in from the read 
statement labelled ONT as long as the record counter is not zero. Also 
notice the two statements before the label ONT. These two statements ef- 
fectively skip the first twenty five lines of each readout (when REC=1) 
for reasons which are stated later. The program segment from the label 
MINN to the label THRTY is probably the most difficult part of the entire 
program. As end of file and parity problems x*ere encountered, it was ne- 
cessary to revise this section continually. The logic shown is a combi- 
nation of that which was needed for routine execution of the program and 
that which was needed to cover any abnormalities which arose. Initially, 
the final twenty five line record of each readout was partially filled 
with anywhere from five to twenty five lines of data (in groups of five 
lines) while the rest of the lines in the record were filled with zeros 
(it was possible for this last record to be completely filled as it turned 
out). This incomplete record was being used initially to indicate the end 
of a readout but the program was eventually modified in order not to rely on 
this. As a result of the change though, the last partially filled record 
of each readout was not written onto the output tape. The third line before 
the label OAHU checks the satellite height. If zero, this would signal the 
end of a readout. The variable X would be set to one and this variable 
along with E would be passed back into the main program where it was tested. 
The segment from the label OAHU to the label FIVE sets up the ORBIT array, 
checking the condition of each line of data and replacing bad lines with 
9999. The loop from label FIVE to label SIX sets up the black and white 
plates and cones data in the SENSOR array, and subtracts 1024 from any count 
values which are greater than 1024, The loop from the label MAUI to the 
label SSP performs the same operations for the mount counts and the refe- 
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rence resistor count. The next program segment from label SSP to label 
OKLA was a later addition to the program. Errors in time were detected 
and they showed up in the zenith angles of the day night transition zones. 
These errors were found to cause as much as a nine degree discrepancy in 
the subpoint latitude. This made it necessary for every readout to be 
checked for this problem. As stated, Record one of each readout is always 
skipped since experience showed that data in this record were often ill 
conditioned and were eliminated anyway. It was decided to allow the mecha- 
nics of the data system settle somewhat, before this time testing would 
be performed. As a result. Records three thru nine were tested for this 
possible error in time. These seven records yield 175 lines of data and 
this corresponds to approximately 85 minutes of orbit time. In checking 
for this error, particular attention is paid to the transition zones and 
one of the tested records will most certainly contain such a zone. Speci- 
fically, the test is as follows. If the zenith angle (called D in this 
segment) is less than 119 degrees then the satellite should be in the sun. 
White plate counts should then be greater than 350. If these counts are 
less than 350 a time error is indicated along with the record number where 
it was detected. This indication appears in the printed output log. The 
file is immediately closed on the output tape and the rest of the readout 
is skipped. If the zenith angle is greater than 132 degrees then the satel- 
lite should be in darkness and again the white plate counts are tested. 

Here, if these count values are greater than 350 a time error is recorded 
as before. Since this test had to be added to the program, partial readouts 
now occur on the output tapes. These partial readouts can be skipped by 
noting the skipped orbit numbers in the printed log and then bypassing these 
orbits during further data processing. In a very general way, this is all 


53 


that the input procedure does. The remaining segments of the procedure 
are composed of the end of file and parity action labels. There is only 
one remaining test to be discussed. Tape number 16 of the ESSA 7 series 
contained a readout consisting of approximately one 25 line record. A 
test was devised which would determine approximately how far along into 
the data the program had prog7:essed. It was necessary to determine 
rather roughly the orbit number of the last readout on each tape. This can 
be done rather easily if the satellite launch date as well as the beginning 
and end date of each tape are known. For the ESSA 7 tapes these dates 
were written on each of the input tape boxes. 
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A 


25 lines processed and 
written onto tape* 
(Part A) 


2. v 

next 25 lines 
are read in. 



Steps one thru four are repeated continuously. 


Figure Al. Two block approach to data processing. 
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InTEjEh array i-T TS 1 1 ilbOiJ 

Add AY TEddt i : it>o i : 7 ] »TH[o;7 .Oi 1UOO J »Ti ruoi7*0s5o0 J J 

oef I mt rcri»i?)= Tdr T?»T1 CT11»T12)- 

f H t E G E H Alii 

INTEGER y»Z*lF£»# 


I m T E i» t d I»J»iWL* w #PJ 
c OR j^t thru iso no begin 

HF aDiF ILES*FhT 1 » OTTSlUc f OR J*l THRU 7 OU TEMPL!»JJ)i 
IF I FUL 1 I’HtN GO ttj M. ELSE 
AB* LT TStl-J J + to? 

IF ClTbEU NMji Ad THtN lu NJ ORE ELSE 
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hl: e n . ) ; 

for 1*1 truuj r / un 

K t U I ■') 

Fun L«-l M < J l') OH 
UEtflN 

m <■ i u ^ [ + cL-n; 

LJf-L f 

Vv . 1 < (’J-l.JW 

FuK J<*l FriKU ( DO 
LEU IN 

if lol o.o men an tu eng el.se 

IF n ifn H)*K) ThiM ii.i T [] SCOT ELSE 

IV It P'il + i»JJ t )L 0 . u THEN GO VO tNG ELaE 

T | H . U > <• T t M P L i * J ] h ■■ x ( TEMPI 1 + 1 . Jj"TEHPL [ #0 J ) i 

Gij TO l J u * 

SCUT * ri* M-] 000; 

if eol o.o then an in eng el&e 

HCk.jT* TENP l T . j J V V x C TEMF U + l . J ] -TEMP U . J 1'J I 
f n a j enu; 

Loot 

LOO# 

Frio 1*1 THUU 1 !)u H f a 1 w 
LUCK L Tp t 1 ,* J ) t 
LUCK IT IK L n *i)i 
ENOi 

CLOSEIFIlEoJ; 
o*'jyiJv,oiifiE*.«y/ii 
SIXTH! KECt-i: 

THIRDS MV!f 
F InST i eoEC + i j 

INPUT ( ‘ L- . aUkMT , AtjFoS.JII . E > X 1 ; 

IF A FiL ^ IMLM at) Til b I X 1 rl ‘■'LSt 

IF t ENL <L THLn SU Hj UHE 'ELSE 

i ► kEc t jl o Then an m Kinar flSl 

SKC'Jnu * Kfc O NEC + 1! 

1 ■ i F I J T ( r< F L » J 0 K o i T . n S F 0 S Li H » F . \ ) ; 

!>■ X F pl ? THt'N 0.1 TO K1XTH Fl.SE 
iV E EaL d ThEh Oil TO URL LL$ r 
P h K T A : 

FuK J*o ihku 2 3 nn 
h e a i n 

L A ri t L .i • C HS 1 K * S C R * P I T T ► P H I L A • P UR T * 


i V aSEhSURL J. 1 J LOL 0.0 l’ HEN SO TU IiJWa ElSl 
IF ASENSOftt j. 2 J Fat oyyy.o rnEN Go TU PrilLA ELSE 
IF AUkuITl J* 4J LSS U1.U CHFN GU TU PlTT ELSE 
It A S t < • S U i< L J + 1 . ? 1 F J L 0999.0 THEN GU TU CHSTR ELSE 
u t u T a f * aoensup t j+ 1 . ? j - a Sen surf I J . a J J! 

OELANG* AiJrttJi r ( J+l , A l-AUKO t T 1 ’ JM ] ; 

!>ElTAG<- AdSl ASE>SUPt J. 2T “ASFNSUtU J> 1 J ) 5 
IF M f oL 1 I'HtN ail TU SCK else 
It- UtLfAf ulK 100. 0 THEN H<-1 ELSE bu TU WHi 
A SENS HI I <■ ASFmSiJKI J»2J + 10.34.0; 

an in PH I la; 

/tiil IF Uf.LTafj LSS loo . 0 I HEN H«-l EL^E GU TU PITTS 
nil I U »- H I L A f 

CmSTH: F i. \ t\«-j JhKtJ A DU ASEobmU J»K J4- 9999.0; 

G j TO m hIla; 

SCrti IV SluutDLLAOsil arn -1 IHFN GU FU PhILA ELSE 
lV ULLTAV LSS -100. :1 THEN Mt-0 ELSE GO JIT PUrT; 

an r ti ila; 

pijin; if AStNbuia j. i i .a n< 350.0 then m<vo fl^e arj to Phila; 

jll Ty FrllL*? 
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»-» L IT : Abi . Nbiji* i Jf !> i t Abf''JSUtU J * 2 1 + 1024 . 0 ; 

,J * 1 ! L A * Jf i Jt 1 t : 

F.nu* 

It- Aih-jiaK: J4 »l , p ; ,: MEN Oil fu NY ELSE 

l\ Ail [ ft .i-tf 'I i L l ii llb.o f H F i j J T U mi rt ELSE 

It- tj j f ■ . b i j t i k ■ * / i? J t.-M. v99‘?,o TH' \ *0 T'J huh ELSE 

ty t. L I <* f- «• ; * a ■_ f-| a i i P 1 1 » * 1 J ■" m 'a t. 'o 0 R i 24*2]? 

Ot.Lt - n< Mis '.a [ I ; J , «. i -4ijN.il f L?.t * 4 1 J 
i 1 1 1. t - b t ' b > ! • f * A * 2 1 “ A S *■ n b ' J N L 2 *4 * 1 J > ? 

il- «■' fJL I 1 ml r. u fij bCH El. :.E 

IF t,i L i i 1 1< U(i» i f rl*N bubE 'ilJ Tit U(iH? 

A S l ' < b •. t "* l 2 *’ * 2 ) A S f 1 5 1 ) H L 2 4 t 2 1 + i 0 2 4 > 0 i 

(»ij li> 'vY? 

H'jhi It- ut-LlA>i Lib 1 > 94 „U iTIcW MM Lube' oJ T U ALbJ 
(ill Hi nY? 

tHjKt : Fup l\tu I H *i l j 4 DU ASENbiH C 2-4 * K ] <• 9999*0? 

G U Hi m i ? 

SCH! H 5 t ON L OlL AMU ) ijTR -] i HFN ulJ TO NY ELSE 
I> OLLlAF lSS “100*1 THEN >Hn ELb£ GiJ TO ALO? 

on i u ny; 

A L U i 1 1- AS £ -JSuH L ?4 , I | OTN ibO.O THEN »-l *■ 0 FLSE GO TO NY? 
uO r u t4 y ; 

aLO: aS£nS‘]IH d A * 2 J t a aFMSUH l 24 * ? ] + 1024*0? 

NY i ASENbriKi24»bJt 1 1 ; 

T U H :,*0 T.iHU 4 o t 
t‘ r. t» 1 1-! 

I. A M £ L H li X * L L M ? 

Ft £ x d{ 

lr MbLi' ilim Ft t * 0 I > OL 0*0- T H E i\j uO fu MJ Y ELSE 
K 00 lt t > f :1 H ‘J 4 OD 
H £ 'j A > . 

It- AbL NbUttlF + i *01 1 JL 4994,0 f He N OU TO ELM ELbL 
1 1- *.w.Lwbttr: LF + i * :i 1 . Sb 4 0 0.0 THE>x UU Til tsKX ELbE 
G lj III U L ,v l r 

SAX.: Abt ' -bi| 4 ;.F + i»uit A-iFlUbDKtF + I >i)l + 102 / 4 . 0 ? 
r ui-l ! ENl • 
t- U 0 i 


pa Y : i- uH i. to r 1 1 1* u in o 

ritlli 1 O' 

L A ij I. L bLMfH.M'nU L»4 . i 


t Jli Jf I Tm <U 4 t )'| 
f E M 1 .. 

C t bt A b P 4 HiJ H L I r J 1 I 



IF 

L 1 

S L‘. ii l. 

0 IKES 


IK G 

1 b 

LUL 

ft I M K N i 


If 

il S jTii 

1 U U 0 , 0 


AoL 

-ib 

: J 4 L 1 * 

JJt TiC'l 


GO 

1 P 

H 4 1 I 


S'.H 

» g i 

j t 

Abb'ibUK L ) * J I 


H K 

IS 

r o Tri 

aoli THF ; 

n u : 

A C5u 


TUli 

JJt Tl[c 

ukT 

: Eg 

u ; 



Tl.L 

; l 

■•iii 

A 

t 


TLX 

: K 

o u 

t o 

TupO 4 ; 


»u ** 

i n 




L«UK 

V. 

t. A * ;JH 

A fl 0 ? 


r-iCb ,944? 
il TO TLX FLbL 
I TU ILL FL St: 

■■HEN GO TO bLO ELSE 
5 » J J • 

“ 1000 . 0 : 

CT ? Fl SE LiU TU AG ? 
r b» J ! ? 
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FQK 1*3 TrtKU 4 0 1 ] 

BEGIN 

CTo<* Abe HbUhLR+I 1 1 ] TIriEb ,994; 

1 F LTS Eol U f H L 4 >,U ID NU El.bE 

U' cib t il H then u>i ru n h else 

T F CTb **Tit loOo.o ThEh GlJ T:1 LA ELSE 
AoEh&tM tH + Wo J* TLO f 3 j* I + 3 J i 
GU Tu .M ;1 » 

(Ai LJS* AbENoutU F + I * 0 1 - 1 uOO . J t 

IF C i S 1 F < Wi H'M CFb^SoJ EL.bE GU TO HUi 
Wut AoEfoO^U-' + i »u 1 * THCTa,I + j|; 

NHi F*|0; 

End; 

NJM: WKlTttur » <r»ftET fF[]t< [4-0 THRU ?4 do LfQ|? J*0 THRU 9 DO AGR9 iT Cl > 0 J * 
F'UR K*0 THRU «j DU- AStNSOfU 1 »K J i 3 1 
IF X wE'l 1 I'ME.M Of) TO STE ELSE 

" CEu5'EnnT*j) * 

uu tu sixth; 

5 re; kETC <- 4tc + i; 

INHUT (REC. AUHdlT. AaEi'iSUKfE’iX) ! 

ir X Eul. 3 THEN GU TO SIXTH HL.St. 

IF E E‘JL 2 fiitw GO TO UHE ELSE 

’PAN'TBT" ' ’ 

F UK J*U THRU 33 00 
BEGIN 

L AdFlL UC pSrt »UALS » i\UST*HS TN » Sl)i 

IF 3 SENS )»U U p t J ■_ 0 L 9.0 THEN GO fO Ubvt ELSE 
IF "bSENdUKt J*31 FOl 9999,0 THEN Go TO SU ELSE 
IF U 1 1 H H I riJuAJ LSS llb.O THEN GU TO hstn else 
IF b b fc, 0 S i j rt l J + 1 ,31 HOI. 9999,0 THEN GU TO DAlS ELSE 
UEu T aF *f UbLHbUR C J + 1 * 3 i **0 bENSUR L J , 3 i i 
OELANG* BUl<unU + l,«] “ HOHOITtO *4 I ; 

ULuT AG<- Abb t ribENS JH l J * 2 J “dSENSOH [ J . 1 1 ) J 

IF u Ehl i Then go to aust hlse 

IF UELTaF G f H 100.0 THEN '1*1 ELbE uU TO ba* 
bbENb JHL J»3 J* BSFNGnR £ 0 1> 2 J + 1034.01 
UO TO o'); 

sa; u' delta* uba 100 . o then m<*i else gu tu hstn; 

gu tu su; 

TdALSI " T OK K*0 ThHU a 1)0 bSENSURt J»M* 9999. Of 
Go T 0 SO i 

rtUb r: IF b IGwdJELAiJG ) GTk “1 THEN uO TO SO ELSE 

iF UELTaF Lbb “loo.o the I H*0 ELSE Go ru uc; 
so ro so; 

HCi IF uSEnSIKI J. i ] GTN 3S0.0 THEN 0*0 ELSE uO TO SD; 

GfJ Tu bu; 

H S T N • HbLNbOH L J»3 J * HbFNSUK L J . ?. J + 1034.0; 

SUt UbENSOHt J.5i*H! 

enu; 

IF hbEuSiJK[34« ? 1 LOL 9999. U THEN GU T(1 SAC ELSE 
jh hUKri i T L 34 * -4 1 L-sfi 11^.0 THEN uO TlJ CAL ELSE 
U AbENbOKLu.d] F<n 9999, u THEN GO TU AHN ELSE 
1 ‘ELTaF^ AbENbuR C 0 » ? ] ** HSLNSOR L 24 . 3 J \ 

' ULLaNG* AUf'b I T | 0 ,4 l-HUKbl T [?4»4 J ; 

UEL1 AG* «US l HoE NSUR l 34p3]“USFNSUrtL34» 1 ] ) ; 

If h LOL. i then GO If) NACij FLbE 
IF UELThI blR Tuu.O ThtO ilt-1 ELbE GU TU Abi; 

Hbt\bi)Hi 34.3 I* 0 SLObllK C 34 . 3 1 + 1034 , Uf 
G IJ 10 b,HC; 

Abl; lF OtLiAu Lbo luu.n ImF.N M4-1 else GO TU CALT 
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ten IQ SmC; 

ATtfi r TUIT KfO THRU 4 iJ'i FiSENSljftl. P4*K 1* 9999*0? 
teij TO SwC;_ 

WAtetj!* ' IF SiteN'rotLAOtiJ teTH “1 ThFN UQ TtJ SAC ELSE 
IF DLLTmS LbS “100.0 THEN H*0 ELSE teU TU SASK# 

GO'Tu SaC; 

SASK« IF tiBENSUKl24,l3 teTR 3BU.0 THEN Ht-U ELSE teQ Tl» bAC# 

teU To Sau; - 

cal: usthsnnt;24>2.i + n£»knsukc?«*? i + loH4.o: 

SAC : BStNblJKi 24*5 mu 
liSVt : FUR L<’U THHu 4 DU 
B £ G IN 

L AbEL SSpUEhV# 

F* 5 * L! 

TF bStNSHKtH + 1 ,01" FQL 0*0 THEN GO TO aYS ELSE 

Fur *<-2 thru a no 
begin 

IF USENteORCH + UOl FQL 9y9V.O THEN «fi TO OEhV ELSE 
IF bSEn'S'oht'P+lJoi t.ss aug.o then ou Tn ss else 

teUTUUENVI 

S5T BSF.NSURCP + IpOH bVfNSURCP + I ,0 J + 1 G R 4 a 0 # 

oenv: end? 

E NOT" 

AYS: TqfT'Uo Thru ?n no 
HEtilN 

LABEL VA * WaSHVINO • Bf? J 

Fqr j*i thru a nn 
BEGIN 

CTB F HS Ef ■ S'tl K L I V J J TIMES *994! 

IF CTS UiL 0 Then fiO Ttl AHI'Z ELSE 
IF’ CTS EuL H THEN ten TO 1NU Ft SE 
IF CTS teTR 1000*0 THEN teQ Ttl VA ELSE 
BSEnSURU * J 1* T t C T S * J ] * 

GU Ul naSh? 

va: CIS* uSEnSuk [I *U 1 - 10U0.0? 

IK CTS ill R 500 THEN CT3<-500 n. SE Go 10 BK# 

HR: B SEN SON 1 1 » J 3 <• TlfC*fSpJ3; 

wash: END ? 

“ TND : ’END; 

ARIZs FOR Kt-O THfiu a lift 

BEGIN 

LABEL ORL * FLA * 1 Tn; 

R 5 x n ? 

For i <rv thru a on 

BttelN 

CTS* hstf SURIP+T >0 J TIMES .994? 

IF “crs fql o them on to colo else 

IF CTS UL H ThFN teO TU FLa FlSE 
IF CTS i»Ti\ 1000.0 ‘l HEN teO Td ORL ELSE 
bSEN&UKtP + I/'OJ< ru:is»l+3 j; 

Git TO FLa# 

ORL: ' CTb* tiSENSuRLP+t *01-100 O.Oi 
" IF CTS teTP 500 THFN CTS*50U ELSE Go TO TTN# 

TIN! b St NS OR [R+I *01* mCTSfrU3l; 

FI a: 1 19 D i 

END? 

CULCJS RK IT £. i □ T > * *rt£Y * } OR 1*0 THRU 24 00 [FOR J*0 THRU 6 DO BURBITUfJJ 
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FIjK k*u IhHii r . UU hSLNbdrtl 1 »K ] J ) » 
|F X NFO 1 Thtiv un TIT srcodu F'lSL 
C L U £>£. ( U I » * ) i 
$t : 

(>u in six 1 hi 

OHE : 

1 ‘LIjSl i u f * ) » 

IF X Em. 3 THFI. C U ) S f ( 0 T * * ) FLSF. 

lnu ; 


fc> u . 


KKlihb I'LTEtltl) - 0. CUMPTLATUlN 7 I ML = ib SLCUNLtS. 

JL_UV > TOTAL StbMli.T SlZt = 1 1 7 1 WOHOS* DISK SIZE = ol SLGoJ NU. PbM. btBS = 31 

UKE STUPAGL KfuiliKtU = t>8 BP KURDS. 

? 

U> i 1. 1 Y KFXIFY hEuUIKlp = 0 wtJPOS. 

AHD-lMAGLS HKJCL'SSlU = bPO. 



UUOOO/L ILL 731 JaYLlSfcP = HH:U44iU)MHILL JL'h ArEU/RH3 1 44 XmLGOL 

r i 

OJLOO/LINL 73 1 3oVOSt ^ = l'M3 14<4 S LOMP ILL oLmAmEu/Ph J 1 44 X.4LGQL 


□ 



. KUHl 1 * 
• MJHf I # 





origin m.W£«® 
of ?OOft 
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